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A convincing detection of primordial non-Gaussianity in the local form of the bispectrum, whose 
amplitude is given by the /nl parameter, offers a powerful test of inflation. In this paper we calculate 
the modification of two-point cross-correlation statistics of weak lensing - galaxy-galaxy lensing 
and galaxy-Cosmic Microwave Background (CMB) cross-correlation - due to /nl. We derive and 
calculate the covariance matrix of galaxy-galaxy lensing including cosmic variance terms. We focus 
on large scales (I < 100) for which the shape noise of the shear measurement becomes irrelevant and 
cosmic variance dominates the error budget. For a modest degree of non-Gaussianity, /nl = ±50, 
modifications of the galaxy-galaxy lensing signal at the 10% level are seen on scales R ~ 300 Mpc, 
and grow rapidly toward larger scales as oc R 2 . We also see a clear signature of the baryonic acoustic 
oscillation feature in the matter power spectrum at ~ 150 Mpc, which can be measured by next- 
generation lensing experiments. In addition we can probe the local-form primordial non-Gaussianity 
in the galaxy-CMB lensing signal by correlating the lensing potential reconstructed from CMB with 
high-z galaxies. For example, for /nl = ±50, we find that the galaxy-CMB lensing cross power 
spectrum is modified by ~ 10% at I ~ 40, and by a factor of two at I ~ 10, for a population of 
galaxies at z — 2 with a bias of 2. The effect is greater for more highly biased populations at larger 
z; thus, high-z galaxy surveys cross-correlated with CMB offer a yet another probe of primordial 
non- Gaussianity. 

PACS numbers: 98.62.Sb; 98.65.-r; 98.80.-k 
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I. INTRODUCTION 

Why study non-Gaussianity? For many years it was 
recognized that the simple inflationary models based 
upon a single slowly-rolling scalar field would predict 
nearly Gaussian primordial fluctuations. In particular, 
when we parametrize the magnitude of non-Gaussianity 
in the primordial curvature perturbations £, which gives 
the observed temperature anisotropy in the Cosmic 
Microwave Background (CMB) in the Sachs-Wolfe 
limit as AT/T = — £/5, using the so-called non-linear 
parameter /nl Q as CM = Cl( x ) ± (3/nl/5)C£(x), 
then the bispectrum of £ is given by 1 (fci, ki, k^) = 
(6/nl/5) [P ( (k 1 )P c (k 2 ) + (2 cyclic terms)], where 
P((k) oc fc™ 3-4 is the power spectrum of £ and n s 
is the tilt of the power spectrum, constrained as 
n s = 0.960 ± 0.013 by the WMAP 5-year data @. This 
form of the bispectrum has the maximum signal in the 
so-called squeezed triangle for which k$ <C &2 ~ ki Q. 
In this limit we obtain 
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B c (k u ki,ka ->()) = — /N L P ( ( fc i)^c(fc3) 



(1) 



'Electronic address: djeong@astro.as.utexas.edu 

1 Definition of the bispectrum in terms of Fourier coefficients of £ 
is <C kl Ck 2 Ck 3 > = (27r) 3 <5(ki+k 2 +k 3 )B c (fci,fc2,fc3)- Throughout 
this paper we shall order ki such that k$ < k2 < ki. 



The earlier calculations showed that /nl from single-field 
slow-roll inflation would be of order the slow-roll parame- 
ter, e ~ 10 -2 0,[E[1]. However, it is not until recent that 
it is finally realized that the coefficient of PQ(k 1 )P ( ^(k 3 ) 
from the simplest single-field slow-roll inflation with the 
canonical kinetic term in the squeezed limit is given pre- 
cisely by 0,1 

B ( (k u k u k 3 - 0) = (1 -n.)P c (fci)P c (fc 3 ). (2) 

Comparing this result with the form predicted by the /nl 
model, one obtains /nl = (5/12)(l — n s ). 

Perhaps, the most important theoretical discovery re- 
garding primordial non-Gaussianity from inflation over 
the last few years is that, not only models with the canon- 
ical kinetic term, but all single-inflation models predict 
the bispectrum in the squeezed limit given by Eq. ^ty, 
regardless of the form of potential, kinetic term, slow- 
roll, or initial vacuum state [1, [l(| [ll|, OjJ . Therefore, 
the prediction from all single-field inflation models is 
/nl = (5/12)(1-ti s ) = 0.017 for n s = 0.96. A convincing 
detection of /nl well above this level is a breakthrough 
in our understanding of the physics of very early universe 
[IE [3. The current limit from the WMAP 5-year data 
is / NL = 38 ± 21 (68% CL) QJ]. 

There are many ways of measuring /nl- The most 
popular method has been the bispectrum of CMB [l|, 
lm [l7l . [H, [3 (also see [2(| for a pedagogical review). 
The other methods include the trispectrum of CMB [21], 
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and the 



|22j, the bispectrum of galaxies [2J, |2J, |25|, |2 
abundance of galaxies and clusters of galaxies [2 

m. 

Recently, analytical [111, [H, [H, [H, H|| and numerical 
[3ll 13(1 [33, HI] studies of the effects of primordial non- 
Gaussianity on the power spectrum of dark matter ha- 
los, Ph(k), have revealed an unexpected signature of pri- 
mordial non-Gaussianity in the form of a scale- dependent 
galaxy bias, i.e., Ph(k) = b\P m (k) — ► [61 + Ab(k)] 2 P m (k), 
where P m (k) is the power spectrum of matter density 
fluctuations, and 



Ab(k) 



3(6i - 1)/niA 



D{z)k 2 T{k) 



(3) 



Here, D(z) and T{k) are the growth rate and the trans- 
fer function for linear matter density fluctuations, respec- 
tively, and 8 C — 1.68 is the threshold linear density con- 
trast for a spherical collapse of an overdensity region. 
The k 2 factor in the denominator of Ab(k) shows that 
this effect is important only on very large scales. Highly 
biased tracers are more sensitive to Jnl- 



II. HALO-MASS CORRELATION FROM 
GALAXY-GALAXY LENSING 
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FIG. 1: Coordinate system and 71 and 72. The shear along e.\ 
has 71 > and 72 = 0, whereas the shear along e2 has 71 < 
and 72 = 0. The shear along e.\ + has 71 = and 72 > 0, 
whereas The shear along ei — e.i has 71 = and 72 < 0. 



A. Formula 

The scale-dependent bias was theoretically discovered 
when the authors of [3l[ studied the form of the cross- 
correlation power spectrum between the dark matter 
halos and the underlying matter density fluctuations, 
Ph m {k) = [h + Ab(k)]P m (k). We can observe P hm {k) 
by cross-correlating the locations of galaxies or clusters 
of galaxies with the matter density fluctuations traced by 
the weak gravitational lensing (see [39[ for a review) . 

One efficient way of measuring P/ lm (fc) is to use the 
so-called galaxy-galaxy lensing technique [S(| [H, |H, 
44, 45, 46]: choose one lens galaxy at a redshift z^, and 
measure the mean of tangential shears in images of lensed 
(source or background) galaxies around the chosen cen- 
tral lensing galaxy as a function of radii from that cen- 
tral galaxy. Finally, average those mean tangential shears 
over all lensing galaxies at the same redshift, zl- 

We begin with the definition of the tangential shear, 
7i, on the flat sky 2 

It (0) = -7i W cos(20) - 72 (0) sin(2<£) , (4) 

where 6 = (0 cos <fr, 6 sin cj)) , and 71 and 72 are compo- 
nents of the shear field. 3 The coordinate system and 



2 For an all-sky analysis, this relation needs to be replaced with 
the exact relation using the spin-2 harmonics [47|| . 

3 As the shear has two independent components, we are ignoring 
another linear combination of 71 and 72 by only focusing on the 



the meaning of 71 and 72 are explained in Fig. [T] For 
purely tangential shears shown in Fig. [lj 74 is always 
positive. This property allows us to average 74 over the 
ring around the origin to estimate the mean tangential 
shear, j t : 



1M = 




On the flat sky, 71 and 72 are related to the projected 
mass density fluctuation in Fourier space, ac(1), as 

li(0) = J -^ K (l)cos(2^) e il e , (6) 

12(9) = J ^L K (l)sin(2^)e il - e , (7) 

where ip is the angle between 1 and e±, i.e., 1 = 
(I cos ip, I sin ip). Using Eqs. ([6]) and ([7]) in Eq. (TJ|, we 
write the tangential shear in terms of k(1) as 

lt(0) = J 7^2 «00 cos[2(0 - # Wcos (^). (8) 



tangential shear. In particular, on large scales there is informa- 
tion in the other component of the shear, and thus the full anal- 
ysis including both shear components (not just tangential one) 
yields a modest (smaller than a factor of \/2) improvement in the 
signal-to-noise ratio. Moreover, using magnification (in addition 
to shears), which is proportional to the convergence field k, can 
also yield a modest improvement. 
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The mean tangential shear (Eq. ([5])) is then given by 
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k(I)J 2 (W). 



cos[2(^-p)]e <WcGB (* - *'> 



(9) 



Here, we have used the identity 
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i(mip- 



(10) 



with m = 2, ip = (j> — cp — 7r/2, a = <£> + 7r/2, and 
dilisia(2'ip)e ixcos ^ = 0. 

The ensemble average of the mean tangential shear 
vanishes, i.e., (j t ) = 0, as («) = 0. This simply means 
that the average of the mean tangential shears, measured 
with respect to random points on the sky, vanishes. We 
obtain non-zero values when we average the mean tan- 
gential shears measured with respect to the locations of 
halos (galaxies or clusters of galaxies). This quantity, 
called the galaxy-galaxy lensing or cluster-galaxy lensing, 
can be used to measure the halo-mass cross correlation. 

While clusters of galaxies may be identified directly 
with dark matter halos of a given mass, how are galaxies 
related to halos? Some galaxies ("field galaxies") may 
also be identified directly with dark matter halos; how- 
ever, galaxies residing within groups or clusters of galax- 
ies should be identified with subhalos moving in a big- 
ger dark matter halo. For such subhalos our argument 
given below may not be immediately used. However, it 
is observationally feasible to identify the central galaxies 
in groups or clusters of galaxies and measure the mean 
tangential shear around them. A number of studies of 
Luminous Red Galaxies (LRGs) extracted from the the 
Sloan Digital Sky Survey (SDSS) have shown that these 
are typical central galaxies in galaxy groups 0, |H, [4t| . 
Scalings such as the mass-luminosity scaling imply that 
LRGs provide a useful proxy for the halos within which 
they reside. We will assume in this study that such trac- 
ers will enable the halo-shear cross-correlation to be mea- 
sured. There are some caveats such as bimodal mass dis- 
tributions in galaxy groups [49[ and the extrapolation to 
higher redshift, but we will leave a detailed exploration 
to real galaxy tracers for later work. 

The ensemble average of the mean tangential shears 
relative to the locations of halos at a given redshift zl, 
denoted as (7^ )(6,zl), is related to the angular cross- 
correlation power spectrum of halos and k, C/" 1 , as [5(| 



l -^C^{z L )J 2 {W). (11) 



We give the derivation of this result in Appendix [AJ 

With the lens redshift zl known (from spectroscopic 
observations), we can calculate the comoving radius, R, 
corresponding to the angular separation on the sky, as 
R = 6<Ia(0)Zl) where cIa(0;zl) is the comoving angular 
diameter distance from z = to z = zl- Using Limber's 
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FIG. 2: Critical surface density, Ti c (zl; zs), as a function of 
the source redshift, zs, for various lens redshifts that roughly 
correspond to the Two Degree Field Galaxy Redshift Survey 
(2dFGRS; zl = 0.1, solid), the main sample of the Sloan 
Digital Sky Survey (SDSS; zl ~ 0.2, dotted), the Luminous 
Red Galaxies (LRGs) of SDSS (z L = 0.3, dashed), and the 
Large Synoptic Survey Telescope (LSST; zl = 0.5 and 0.8, 
dot-dashed and triple-dot-dashed, respectively). 



approximation [5l], [52j on the flat sky relating C/ 1 " to 
Phm(k), 4 we can write Eq. (flTj) as [5(| 

(^)(R,z L ) = [ ^P hm (k,z L )J 2 (kR). (12) 

Here, po is the mean comoving mass density of the uni- 
verse, and Ti c {zl) is the so-called critical surface density: 

^c 1 ^) = ^f-{l+ZL)d A (0;z L ) 



d A (z L ;z s ) 
dz s p (zs)-j —Tf. s", 

where p(z$) is the redshift distribution of sources normal- 
ized to unity, J dzp(z) — 1, and oU(0;z) and dA^z; zg) 
are the comoving angular diameter distances out to z 
and between z and zs, respectively. The numerical value 
of 4ttG/c 2 is 6.01 x 10~ 19 Mpc/M Q , and 4ttG> /c 2 is 
1.67 x 10- 7 (n m h 2 ) Mpc~ 2 . 
Eq. (fT2| is often written as 



(H)(R,zl) 



AE(i?, z L ) 



(14) 



£ c (zl) 

To simplify the analysis, let us define the "effective source 
redshift" of a given survey from the following equation: 



dAjzL', Zs,ef{) 

dA(0;z s , c s) 



dzs ( /A;*s) 
L dzsp[zs) d A (0;z s ) ■ 



(15) 



4 As we arc dealing with correlations on very large angular scales, 
one may worry about the validity of Limber's approximation. In 
Appendix \C\ we give a detailed study of the validity and limita- 
tion of Limber's approximation for the galaxy-galaxy lensing. 
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FIG. 3: The baryonic feature in the matter power spec- 
trum, as seen in the galaxy-galaxy lensing, AE(i?), for sev- 
eral populations of lens galaxies with 6i = 2 at zl — 0.3 
(similar to SDSS LRGs, solid), &i = 2 at z L = 0.5 (higher- 
z LRGs, dotted), b\ = 2 at zl = 0.8 (galaxies that can 
be observed by LSST, dashed), and 6i = 5 at zl — 0.8 
(clusters of galaxies that can be observed by LSST, dot- 
dashed). The vertical line shows the location of the bary- 
onic feature, -Rbao = 106.9 hT 1 Mpc, calculated from the 
"WMAP+BAO+SN ML" parameters in Table 1 of 0. Note 
that we have used the linear matter power spectrum and the 
Gaussian initial condition (/nl = 0) for this calculation. 



Henceforth we shall use zs to denote zs,cff, and write 

v-lf > 4?rG n , \,i fa ,d A (z L ;z s ) 

S c (z L ;zs) = ^-(1 + z L )d A (Q;z L )—— -. (16) 

Fig.rjshows S c for z L = 0.1 (2dFGRS, Two Degree Field 
Galaxy Redshift Survey), 0.2 (SDSS main), 0.3 (SDSS 
LRG), and 0.5 and 0.8 (both LSST, Large Synoptic Sur- 
vey Telescope). The smaller E c is, the larger the observed 
mean tangential shear is. 



B. Results 

We can now calculate the observable, AE(i2, Zl), for 
various values of /nl- We use 

AE(i?, z L ) 

/kdk 
——P m (k,z L )J 2 {kR) 
2tt 

/kdk 
— Ab(k 1 z L )P m {k 1 z L )J 2 (kR), (17) 

where the scale-dependent bias, Ab(k,z), is given by 
Eq. ([3|). As we are interested in large scales, i.e., R > 
10 hr 1 Mpc, we shall use the linear matter spectrum for 

Pm{k). 
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FIG. 4: Imprints of the local-type primordial non- 
Gaussianity in the galaxy-galaxy lensing, AE(i?), for the same 
populations of lens galaxies as in Fig. [3] The solid, dashed, 
and dotted lines show /nl = 0, ±50, and ±100, respectively. 
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FIG. 5: Fractional differences between AE(i?) from non- 
Gaussian initial conditions and the Gaussian initial condition, 
|AE(iZ; /nl)/AE(.R; /nl = 0) — 1|, calculated from the curves 
shown in Fig. [4j The dot-dashed, dashed, and dotted lines 
show /nl = ±10, ±50, and ±100, respectively, while the thin 
solid line shows oc R 2 with an arbitrary normalization. 



Fig. [3] shows, for the Gaussian initial condition (/nl = 
0), AE(R,z L ) from R = 50 to 200 h' 1 Mpc. We have 
chosen the bias parameters and lens redshifts to represent 
the existing data sets as well as the future ones: bi = 2 
at zl = 0.3 (similar to the observed values from SDSS 
LRGs [H, top-left), h = 2 at z L = 0.5 (higher-z LRGs 
[H]], top-right), bi = 2 at zl = 0.8 (galaxies that can 
be observed by LSST, [55|, bottom- left), and b\ = 5 at 
zl = 0.8 (clusters of galaxies that can be observed by 
LSST, bottom-right). While LSST is an imaging survey, 
we assume that we can obtain spectroscopic redshifts of 
some (~ 10 6 ) lens galaxies by follow-up observations. It 
is also straightforward to extend our analysis to lenses 
selected by photometric redshifts. 
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At R r-j 110 hr 1 Mpc we see a clear "shoulder" due 
to the baryonic feature in the linear matter power spec- 
trum (often called Baryon Acoustic Oscillations; BAO). 
The sound horizon at the drag epoch (which is more 
relevant to the matter power spectrum than the pho- 
ton decoupling epoch for the CMB power spectrum) 
calculated from the cosmological model that we use, 
the "WMAP+BAO+SN ML" parameters in Table 1 of 
0, is 106.9 h- 1 Mpc, as shown as the vertical line 
in this figure. The magnitude of AE on this scale is 
~ 0.1 h Mq pc~ 2 . Assuming a range of E c from future 
surveys, E c - 1000 - 4000 h M Q pc" 2 (see Fig. HJ}, this 
value corresponds to the mean tangential shear of order 
2.5 x 10~ 5 to 10~ 4 . Is this observable? 

For comparison, Sheldon et al. measured AS(iZ) ~ 
0.5 h M Q pc~ 2 at R ~ 30 /i -1 Mpc from clusters 
of galaxies in the SDSS main sample. The mean lens 
redshift for these data is zl ~ 0.2, which would give 
E c - 5000 h Mq pc~ 2 (see Fig. [2] for z L = 0.2 and 
zs ~ 0.4); thus, the magnitude of the mean tangential 
shear that they were able to measure is of order 10~ 4 , 
which is only ~ 1 to 4 times larger than the magnitude 
of the signal expected from the BAO. Therefore, detect- 
ing the BAO signature in AE(i?) should be quite feasible 
with the future observations. We shall give a more quan- 
titative discussion on the detectability of BAO from the 
galaxy-galaxy lensing effect in Sec III Dl 

How about /nl? As expected, the effect of /nl is 
enhanced on very large scales, i.e., hundreds of Mpc 
(see Fig. SJ. For / NL = ±50, AE(i?) is modified by 
10-20% at R ~ 300 /i" 1 Mpc (depending on h and 
zl', see Fig. [5]). The modification grows rapidly toward 
larger scales, in proportion to R 2 . On such a large scale 
(R ~ 300 hT x Mpc ), the galaxy-galaxy lensing signal is 
on the order of AE ~ 0.01 h Mq pc -2 , and thus we need 
to measure the mean tangential shear down to the level 
of 7J 1 ~ 2.5 x 10~ 6 to 10" 5 , i.e., 10-40 times smaller than 
the level of sensitivity achieved by the current observa- 
tions. Can we observe such a small shear? 



far as we know this formula has not been derived be- 
fore. Note that we have assumed a single source and 
lens redshift. For multiple source and lens redshifts, the 
covariance matrix needs to be suitably generalized. 

Here, C\ l and C" are the angular power spectra of 
the lens halos (galaxies or cluster of galaxies) and n, re- 
spectively, and nj, and rig are the number densities of 
the lens halos and the lensed (source) galaxies, respec- 
tively. These angular power spectra, Cf K , C/ 1 , Cf, will 
be related to the corresponding three-dimensional power 
spectrum, P(k), in Sec IIII Cl 

In the limit that the cosmic variance is unimportant, 
we recover the usual expression used in the literature: 



mome'))-(it(o))(ine')) = 



S D (0-9') 



'■y D 



N L 2tv9ti<. 



(19) 



where Nl — 47r/ s k y nL is the total number of lens halos 
available in the data. In this limit the errors in different 
radial bins are uncorrelated, and they are simply given 
by the shape noise, cr 7 , reduced by the square-root of the 
number of source galaxies available within each radial bin 
and the total number of lens halos that we can use for 
averaging the mean tangential shear. In particular, at 
each bin with a width A9, we find the variance of 



Varfr^)] 



2Tv9{A9)n s N L ' 



(20) 



in the absence of the cosmic variance. 

When would the cosmic variance become important? 
There is the maximum surface number density of sources, 
«s,max = cr 2 /C; K , above which the shape noise becomes 
irrelevant. This gives the maximum number of sources 
within a given radial bin of a width A9 (<C 9) above 
which the shape noise becomes irrelevant: 



N s<1 



2Tv9{A9)n s , 1 



AO 



(21) 

For 10 = tv (the usual relation between I and 9) and 



C. Covariance matrix of the mean tangential shear cr 7 ~ 0.3 (realistic shape noise), we find 



In order to study the feasibility of measuring the tan- 
gential shear of order 10 -6 , we compute the covariance 
matrix of the mean tangential shears averaged over Nl 
lens galaxies. As derived in Appendix [B] the covariance 
matrix of the mean tangential shear is 

(l1(9)l1(9'))-{^(0)){lt(O r )) 
1 



l ^Lj 2 (l9)J 2 {l9') 
Ztv 




(18) 



N 



S^max 



A9\ 



J 1 2 C?I(2tv\ 



(22) 



At I - 100, 1 2 C?/{2tv) - 10" 5 1(1; thus, we do not 
gain sensitivity any further by having more than, say, 
10 4 galaxies (for A9/9 = 0.1) within a single radial bin. 

Alternatively, one can define the minimum multipole, 
^min, below which the cosmic variance term dominates: 



l2Tvn s l 2 Cf 



2n 



(23) 



For LSST, we expect to have the surface density of 
sources on the order of ns = 30 arcmin~ 2 = 3.5 x 
noise of lens halos, as well as the shape noise <7 7 . As 10 8 sr _1 . For er 7 = 0.3, we find / m i n (LSST) ~ 1.6 x 



This expression includes the cosmic variance, the shot 
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FIG. 6: Same as Fig. [3] but with the expected 1-a uncer- 
tainties for full-sky lens surveys and a single lens redshift. 
Adjacent bins are highly correlated, with the correlation co- 
efficients shown in Fig. [7] The open (filled) boxes show the 
binned uncertainties with (without) the cosmic variance term 
due to the cosmic shear field included. See Eq. (|28l) and (|29[) 
for the formulae giving open and filled boxes, respectively. We 
use the radial bin of size AR = 5 h~ x Mpc. For comparison, 
we also show AE(i?) computed from the smooth power spec- 
trum without the baryonic feature [56[ (dashed lines). Note 
that the uncertainties are calculated for a single lens redshift 
slice, and thus they will go down as we add more lens redshift 
slices. 
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FIG. 7: 
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The cross-correlation-coefficient matrix, 



Cij I \fCiiCjj, where dj is the covariance matrix given in 
Eq. 123, for a radial bin of AR = 5 h^ 1 Mpc. We show r i3 
for the same populations of lens galaxies as shown in Fig. [3] 
and [6] We use the same number of source galaxies and the 
same shape noise as in Fig. [5] The neighboring bins are highly 
correlated for AR < 10 h' 1 Mpc. 



ioV^ 2 cr/( 27r )- At 1 - lo3 > l2 CN( 2lT ) £ io -4 [H; 

thus, at I < 10 3 the cosmic variance term dominates. 

In the limit that the covariance matrix is dominated 
by the cosmic variance terms, we have 

{i h de)i h t {e'))-{i h t {e)){i h t{o')) 

47r/ sky J 2tt 



where n = C/ lK /y CfC/ 1 is the cross-correlation coeffi- 
cient. The variance at a given radial bin is 

Varp#(0)] 

= ^ / l ^[J2(l9)) 2 C?Cni+ rh (24) 

D. Detectability of the mean tangential shear 

In this section, we shall calculate the expected uncer- 
tainties in radially binned measurements of the mean tan- 
gential shear. 

The mean tangential shear averaged within the i-th 
bin, (j t )(6i), i.e., the mean tangential shear averaged 




100 200 300 400 500100 200 300 400 500 
R [Mpc/h] R [Mpc/h] 



FIG. 8: Same as Fig. 21 but with the expected 1-a uncer- 
tainties for full-sky lens surveys and a single lens redshift. 
Adjacent bins are highly correlated. The open (filled) boxes 
show the binned uncertainties with (without) the cosmic vari- 
ance term due to the cosmic shear field included. See Eq. (|28p 
and (|29[) for the formulae giving open and filled boxes, respec- 
tively. We use logarithmic bins with AR — R/10. Note that 
the uncertainties are calculated for a single lens redshift slice, 
and thus they will go down as we add more lens redshift slices. 
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within an annulus between #i !m i n and 6i tXaax , is given by 
2tt 



Aft) 



Odd (>tf)(0) 



(25) 



where = 7r(0f max — 2 min ) is the area of the annulus, 

and 



MlBi) 



2tt 



9d9 J 2 (W), 



(26) 



is the Bessel function averaged within a bin. 

Similarly, the covariance matrix of the binned mean 
tangential shears is given by 

Cy EE $(6i)%m-<n{0i))$V3)) 



1 



47r/ sky J 2ir 

(Q 



3) 

l ^J2mMie 3 ) 



1L 



ns 



.(27) 



ns = 3.5 x 10 8 sr _1 . As the covariance matrix is dom- 
inated by the cosmic variance terms, the size of open 
boxes is insensitive to the exact values of Nl, c 7 , or 
ns- (See Sec. IIII Cl ) First, we calculate the binned un- 
certainties in the region close to the baryonic feature, 
R ~ 110 hr 1 Mpc. In Fig. [6j the open boxes show the 
full uncertainties including the cosmic variance due to 
cosmic shear (Eq. (f2"5)l ). while the filled boxes show the 
uncertainties without the cosmic shear term (Eq. I|29p). 
The latter is clearly negligible compared to the former on 
large scales, R > 50 h^ 1 Mpc. 

Can we distinguish AE(i?) with and without the bary- 
onic feature? Without baryons, we do not see any fea- 
tures in AE(i?); see dashed lines in Fig. [6] which are cal- 
culated from the smooth linear power spectrum without 
the baryonic feature (56|. To see if we can detect this 
feature in AE(i?), we estimate the \ 2 difference between 
AE(i?) with and without the baryonic feature: 



A X 2 = J] (AS, - AE l , nw )C^ 1 (AE J - AE J>W ), 



This matrix contains the full information regarding the 
statistical errors of the binned measurements of the mean 
tangential shear, which includes the cosmic variance er- 
rors due to the cosmic shear (Cf), clustering of lens 
galaxies (C/ 1 ) and their correlations (C7j), the finite 
number density of lenses, and the noise in intrinsic shapes 
of source galaxies. 

The variance at a given radial bin is 



Var[f (0,)] 



47r/ sky J 2tt 



Idl 



[J2(W l )f 



(Cf K ) 2 + 



c\ l + — 

n Lj 



Cf + — m 

ns 



In the analysis of the galaxy-galaxy lensing effects in the 
literature, the cosmic variance due to cosmic shear is usu- 
ally ignored: 



Var[f t (&)] 



1 



K=0 



4vr/ f 



sky 



c l + 



Idl r 

2^ 



Midi 



i 



TIL 



(29) 



This is probably a reasonable approximation for the cur- 
rent measurements at R < 30 h^ 1 Mpc; however, on 
larger scales which will be probed by the next-generation 
lens surveys, the cosmic variance due to cosmic shear 
must be included, as we show in Fig. [6l 

For estimating the expected uncertainties, we assume 
a million lens galaxies with very narrow (delta-function 
like) redshift distribution centered at zl (Nl = 10 6 ) over 
the full sky / s k y = 1. We also assume u 7 = 0.3, and 



where AE^ is the mean tangential shear of z-th bin, AE nw 
is AE without the baryonic feature, and C^ 1 is the in- 
verse of the binned covariance matrix (Eq. |2"7| . Using 
only a single lens redshift slice, we find A% 2 = 0.85 
(z L = 0.3, b = 2), 1.07 (z L = 0.5, b = 2), 1.32 (z L = 0.8, 
6 = 3), and 1.34 [zl — 0.8, 6 = 5). For example, if we add 
up all these measurements at different slices (zl = 0.3, 
0.5 and 0.8), significance of detection of the baryonic fea- 
ture is Ax 2 = 3.2, i.e., 93% C.L. As we expect to have 
many more lens redshift slices from the future lens sur- 
veys, detection and measurement of the baryonic feature 
in AE are quite feasible. For multiple lens slices the gain 
in the signal-to-noise ratio will be approximately \/N\ cns ; 
thus, for 10 lens slices the errors would be a factor of 3 
smaller. At best we can expect ~ 25 slices, which gives 
a factor of 5 reduction in errors. 

What about /nl? We show the expected Tex uncer- 
tainties for the mean tangential shears, AE(i2), on larger 
scales in Fig. [8] For this figure we use logarithmic bins 
with the radial size of AR/R = 0.1. We find that AE(i?) 
on R ~ 250 hr 1 Mpc is detectable, even from a single lens 
redshift slice. This is remarkable; however, the predicted 
uncertainties are too large for us to distinguish between 
/nl = and /nl = 100 using a single lens redshift slice. 
In order to obtain a tight limit on /nl, we would need to 
include many lens redshift slices. 

Note that the uncertainty at a given R is larger for a 
smaller lens redshift. This is because a given R corre- 
sponds to a larger angular size for a lower lens redshift, 
making the cosmic variance contribution greater. 
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FIG. 10: Imprints of the local-type primordial non- 
Gaussianity in the galaxy-convergence cross power spectrum, 
1(1 + 1)C; /(2tt), for for the same populations of lens galax- 
ies as in Fig. [3] The solid, dashed, and dotted lines show 
/nl = 0, ±50, and ±100, respectively. 



FIG. 9: Angular power spectrum of the galaxy-convergence 
cross correlation, C; , at various multipoles as a function of 
the lens redshift, zl, for two effective source redshifts, z s — 1 
(top) and 2 (bottom). We have divided Cj by its maximum 
value. The solid, dotted, dashed, dot-dashed, and triple-dot- 
dashed lines show I — 10, 50, 100, 350, and 1000, respectively. 



III. HARMONIC SPACE APPROACH 



Again using Limber's approximation (whose validity 
and limitation are studied in Appendix[Cj , we find the re- 
lation between the angular cross-correlation power spec- 
trum of the convergence field and the halo density at 
a given lens redshift zl, CJ ik (zl), and the halo-mass 
cross-correlation power spectrum at the same redshift, 
Phm{k 7 z L ), as 



A. Formula 

The mean tangential shear, (7J 1 ) or AS, is currently 
widely used for measuring the halo-shear cross correla- 
tion, as this method is easy to implement and is less 
sensitive to systematic errors. 

In this section, we shall study the effects of /nl on 
the equivalent quantity in harmonic space: the halo- 
convergence cross power spectrum, C, . The mean tan- 
gential shear is related to C/ lK by the 2-dimensional 
Fourier integral given in Eq. pip . 

The convergence field, «(n), is the matter density fluc- 
tuations projected on the sky: 



dzW K (z)S m [d A (0; z)n, z] 



(30) 



where 5 m (r, z) = p m (r, z)/p m (z) - 1, and W K (z) is a lens 
kernel which describes the efficiency of lensing for a given 
redshift distribution of sources, p(zs)- 



W K (z) 



Po 



£ c (z; z s )H(z) 
where the critical density, S c , is defined in Eq. (|16[) 



(31) 



C? k {zl) 



Po 



T, c (zl; zs)d A (0; z L ) 

47rGpo . v 
5 — (1 + z L ) 

x Phm k 



-Ph. 



1 + 1/2 
d A (0;z L ) 



. zl 



d A (z L ;z s ) 



d A (0;z L )d A (0;z s ) 
1 + 1/2 



d A (0;z L ) 



,ZL 



(32) 



Fig. [9] shows C i i IKj (zl) for the Gaussian density field as 
a function of lens redshifts, Zl- The convergence fields 
at low (high) multipoles are better correlated with low-z 
(high-z) galaxies. This is due to the shape of the mat- 
ter power spectrum: on very large scales (i.e., low V), 
the matter power spectrum is given by the initial power 
spectrum, Phm{k) cx k, and thus we get l/d A (0; zl) 
from Phm[k = l/d A (0; zl)]- This gives a larger weight to 
low- 2: galaxies. On smaller scales where Phm{k) oc fc" cff 
with n c ff ~ —3, we get positive powers of d A (Q] Zl) from 
Phm[k = l/d A (0; zl)], which gives a larger weight to high- 
z galaxies. 
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FIG. 11: Fractional differences between C; from non- 
Gaussian initial conditions and the Gaussian initial condition, 
calculated from the curves shown in Fig. [10] These differences 
are equal to \ Ab(l = k/d A , zl)|/&i(z.l). The dashed and dot- 
ted lines show /nl = ±50 and ±100, respectively, while the 
thin solid lines show l~ 2 with an arbitrary normalization. 



B. Result 



FIG. 12: Angular power spectra of the galaxy-galaxy corre- 
lation, C; (thick dotted lines), the galaxy-convergence cross- 
correlation, C; (thick solid lines), and the convergence- 
convergence correlation, Cf (thick dashed lines) for the Gaus- 
sian initial condition (/nl = 0). The four panels show the 
same populations of galaxies and clusters of galaxies as in 
Fig. 1101 We also show the galaxy shot noise, 1/ul (thin dotted 
lines) as well as the source shape noise, a^/ns (thin dashed 
lines), for iV L = 10 6 , cr 1 = 0.3, and n s = 3.5 x 10 s sr" 1 . We 
find 1/ni « Cf and a^/n s < Cf for I < 100. 



We can now calculate Cj for various values of f\ 



NL- 



We use 



— ( 1 + Z L) 



d A (z L ;z s ) 



d A (0; z L )d A (0; zs) 



bi{z L ) + Ab[k= j ~]y 2 , , -L 



d A (0;z L ) : 



x P n 



k = A±±l^ iZL 

<2a(0;zl)' 



(33) 



where the scale-dependent bias, Ab(k,z), is given by 
Eq. ©. 

Figure M shows Cf K (z L ) for / NL = ±50 and ±100 for 
populations of galaxies that we have considered in the 
previous sections. For each lens redshift, we calculate the 
"effective" source redshift by requiring that the angular 
diameter distance to the source redshift is twice as large 
as that to the lens redshift, i.e., d A (0;zs) = 2d A (0;zL)- 
With this requirement, the source redshifts are z s = 0.65, 
1.19, and 2.25 for zl = 0.3, 0.5, and 0.8, respectively. 

Figure [TT] shows the fractional differences between 
non-Gaussian predictions and the Gaussian prediction 
(/nl = 0), which are simply equal to Ab(k, zi)lb\(zh) 
where k = l/d A (0; z£). As expected from the form of the 
scale-dependent bias, the difference grows toward small 
multipoles as roughly 1/Z 2 . While lower redshift popu- 
lations do not show more than 10% difference at I > 10 
for /nl = ±50, a higher-z population of lens galaxies or 
clusters of galaxies at zl — 0.8 show the differences at 
the level of ~ 10% at I ~ 20 and - 30% at I ~ 10. Are 
these effects detectable? 



C. Covariance matrix of the galaxy-convergence 
cross power spectrum 

The covariance matrix of the galaxy-convergence cross- 
correlation power spectrum is given by 



Sw 

(21 + l)/ sky 



(c>n 2 + 




ns 



where Sw is Kronecker's delta symbol showing that the 
angular power spectra at different multipoles are uncor- 
related. Again, C\ l and G* are the angular power spectra 
of the lens halos (galaxies or cluster of galaxies) and k, 
respectively, and Ul and ns are the number densities of 
the lens halos and the lensed (source) galaxies, respec- 
tively. 

We calculate C, K by using Limber's approximation as 



dz 



Pi 



, _ 1+1/2 

K - d A (0;z) 



T? c {z;z s ) H(z)d 2 A (0;z) 



(35) 



However, we cannot use Limber's approximation for Cj 1 
unless one considers lens redshift slices that are broad. 
As we are assuming a thin lens redshift slice throughout 
this paper, we must not use Limber's approximation, but 
evaluate the exact integral relation: 



Ci = - / dkk 2 P g (k,z L )jf [kd A {z L )\ 



(36) 
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FIG. 13: Same as Fig. 1101 with the expected 1-a uncertain- 
ties for full-sky lens surveys and a single lens redshift. Adja- 
cent bins are uncorrelated. The open (filled) boxes show the 
binned uncertainties with (without) the cosmic variance term 
due to the cosmic shear field included. We used Eq. (|34[) for 



the open boxes, and Eq. (J34J with Cf re = = Cf for the filled 
boxes. We use logarithmic bins of Al = 0.23Z. Note that the 
uncertainties are calculated for a single lens redshift slice, and 
thus they will go down as we add more lens redshift slices. 



For the galaxy-convergence power spectra in Fig. [T0l with 
the full sky coverage (/ s ky = 1) and Al/l = 0.23, we find 
l mm = 9.0, 12.1, and 15.7 for z L = 0.3 (z s = 0.65), 0.5 
(1.19), and 0.8 (2.25), respectively. 

Similarly, we can estimate the maximum radius below 
which we can measure the mean tangential shear, AE(i?), 
as 



rcU(0; zl) 



(39) 



For example, with AR/R = Al/l = 0.1, we get i? m ax — 
215, 260, and 300 hr 1 Mpc for z L = 0.3, 0.5 and 0.8, 
respectively. These values do give the radii at which the 
signal-to-noise ratios are roughly unity in Fig. [8j 

Fig. [T3l shows the expected 1-a uncertainties oiCf K for 
several populations of lens galaxies. We find that the cos- 
mic variance completely dominates the uncertainties on 
large scales (low I) where the non-Gaussian effects are the 
largest. Again, while we find that it would be difficult to 
measure /nl from a single lens redshift slice, combining 
many redshift slices should help us measure /nl, espe- 
cially when we can use many slices at moderately high 
rcdshifts. 



where ji is the spherical Bessel function, and P g {k,z) is 
the linear galaxy power spectrum: P g (k) = b\P m (k). 

Fig. [T2] shows the galaxy-galaxy, galaxy-convergence, 
and convergence-convergence angular power spectra for 
Gaussian (/nl = 0) initial conditions. We also show 
the shot noise of the galaxy angular power spectrum, 
1 , and the shape noise of the convergence power spec- 
trum, a^/ns, with the following representative values: 



Nl = iwriL = 10 6 , n s = 3.5 x 10 8 sr" 1 , and ct. 



0.3. 



We find l/n L <C C\ l and a^/ns <C Cf for the multipoles 
that we are interested in, i.e., / < 100, and thus we con- 
clude that the uncertainties are totally dominated by the 
cosmic variance terms. In other words, the size of the 
uncertainties are insensitive to the exact choices of Nl, 
<t 7 , or ns- 

We also find that the va lues of cross correlation coef- 
ficients, n = C^ K I^G\ l Gf, are small (of order 10-20%): 
the maximum values are 0.19, 0.15, and 0.13 for Zl = 0.3, 
0.5, and 0.8, respectively. This implies that one may ig- 
nore the contribution of C7j to the covariance matrix, 
approximating the variance of C/" 4 of a single lens red- 
shift slice for a multipole bin of size Al as: 



Var(C;'" t ) 



{2l + l)Alf sky - 



(37) 



Therefore, we should be able to measure the 
galaxy-convergence cross-power spectrum with 



(38) 



c i K l\ Var(C/ lK ) > 1 when the multipoles satisfy 



l>Jr 



r ly /2(Al/l)U y 



IV. HALO-MASS CORRELATION FROM 
GALAXY-CMB LENSING 

A. Formula 

Instead of using the background galaxies for measuring 
the cosmic shear field due to the intervening mass, one 
can use the CMB as the background light and measure 
the shear field of the CMB lensing due to the intervening 
mass between us and the the photon decoupling epoch at 
z* ~ 1089. See [57) for a review on the CMB lensing. 

The lensing effect makes CMB anisotropics (both tem- 
perature and polarization) non-Gaussian by producing a 
non- vanishing connected four-point function, although it 
does not produce any non- vanishing three-point function. 
One can use this property to reconstruct the lensing po- 
tential field, hence the projected mass-density field be- 
tween us and z*, from the four-point function of CMB 

By cross-correlating the halo over-density field, Sh, at 
some redshift Zl (measured from spectroscopic observa- 
tions) and the k field reconstructed from the CMB lens- 
ing, one can measure the halo-convergence angular power 
spectrum, C/ 1K . 

The angular power spectrum of the galaxy-CMB lens- 
ing cross correlation is merely a special case of the galaxy- 
convergence cross correlation that we have studied in the 
previous section: all we need to do is to set the source 
redshift, zs, to be the redshift of the photon decou- 
pling epoch, ~ 1089, i.e., zg — z*. Note that for 
a flat universe <1a{zl] z *) — rfyi(0;z*) — 6^(0; zjj where 
d A (0;z*) = 9.83 h^ 1 Gpc. 
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FIG. 14: Angular power spectrum of the galaxy-CMB lens- 
ing, C; , at various multipoles as a function of the lens red- 
shift, zl- We have divided C\ by its maximum value. The 
solid, dotted, dashed, dot-dashed, and triple-dot-dashed lines 
show I = 10, 50, 100, 350, and 1000, respectively. 



Figure [H] shows that the CMB lensing at low (high) 
multipoles are better correlated with low-z (high-z) 
galaxies. This is due to the shape of the matter 
power spectrum, as we have explained in the previ- 
ous section. Note that C/ iK of the CMB lensing for 
a given multipole decreases more slowly with z L than 
that of the galaxy lensing due to the geometrical factor 
d A (z L ;z s )/d A (0;z s ). 

Note that CMB and galaxies at Z ?S 1 are correlated 
also via the Integrated Sachs- Wolfe (ISW) effect [p. We 
shall not include this effect in our cross-correlation cal- 
culation for the following reason. We calculate the cross- 
correlation signal between galaxies and the convergence 
field reconstructed from CMB. This reconstruction re- 
lies on the fact that lensed CMB fluctuations have non- 
vanishing connected four-point function. On the other 
hand, the linear ISW effect does not have such a partic- 
ular form of four-point function induced by lensing, and 
thus should not contribute to the reconstructed conver- 
gence field. See [34J for the effects of /nl on the galaxy- 
ISW cross correlation. 



B. Results 

We can now calculate Cj for various values of /nl- 
We use 

nh*( v \ ^Gp Q d A (z L ;z*) 



X P ri 



b 1 (z L ) + Ab[k 
I 



d A (0; ZL)d A {0; z*) 
I 



k = 



d A (Q\ZL) 



d A (0;z L ) 
z L 



ZL 



(40) 



where the scale-dependent bias, A6(fc,z), is given by 
Eq. ©. 

Figure [H] shows Cf K (z L ) for / NL = ±50 and ±100 
for populations of low-z galaxies that we have consid- 
ered in the previous sections: bi = 2 at zl = 0.3 
(similar to SDSS LRGs, top-left), b\ = 2 at z L = 0.5 
(higher- z LRGs, top-right), b\ = 2 at zl = 0.8 (galax- 
ies that can be observed by LSST, bottom- left), and 
b\ = 5 at zl = 0.8 (clusters of galaxies that can be 
observed by LSST, bottom-right). The fractional differ- 
ences between non-Gaussian predictions and the Gaus- 
sian prediction (/nl = 0) are exactly the same as those 
shown in Fig. 1111 in the limit where Limber's approxima- 
tion is valid, the galaxy-convergence power spectrum and 
the galaxy-CMB lensing power spectrum for the same 
lens galaxies differ only by a constant geometrical fac- 
tor of d A [zL\ z*)d A {0; z s )/d A (z L ; z s )d A (0; z*). Inciden- 
tally, for our choice of the source redshifts in the previous 
section, 2d A (zL', z*)/d A (0; z*) = 1.83, 1.73, and 1.60 for 
zl = 0.3, 0.5, and 0.8, respectively. 

Therefore, the galaxy-CMB lensing cross correlation 
would provide a nice cross-check for systematics of the 
galaxy-convergence cross correlation, and vice versa: al- 
ter all, we are measuring the same quantity, Phm(k), 
by two different background sources, high-z galaxies and 
CMB. 

In using high-z galaxies as sources, the galaxy-galaxy 
lensing measurement may be susceptible to systematic 
errors widely discussed in the lensing literature, namely 
shear calibration, coherent point spread function (PSF) 
anisotropy, redshift biases, magnification bias and in- 
trinsic alignments of galaxies. Here we are particu- 
larly concerned with errors that affect galaxy-shear cross- 
correlations by mimicking the angular dependence of the 
signal due to non-zero /nl- Fortunately most systematic 
errors that affect shear-shear correlations do not con- 
tribute to galaxy-shear cross correlations: for instance, 
PSF anisotropy affects background galaxy shapes but not 
foreground galaxy locations [62|. With standard lensing 
data analysis methods, it can be ensured that both the 
shear calibration and PSF do not contribute a scale de- 
pendence to the first order. Biases in the redshift distri- 
butions of lens and source galaxies can similarly lead to 
a mis-estimation of the amplitude of the signal, but not 
its scale dependence. Thus, to the lowest order, the mea- 
surement of /nl via the scale dependence of the galaxy- 
galaxy lensing signal is robust to the leading systematic 
errors in weak lensing. But a detailed study of various 
sources of error is needed given the small signal we are 
seeking. 

Another benefit of using the CMB lensing as a proxy 
for the intervening matter distribution is that we can 
probe the galaxy-matter cross correlation at high red- 
shift to which we cannot reach with the galaxy-galaxy 
lensing method. It is especially useful for probing primor- 
dial non-Gaussianity, as the scale dependent bias signal 
is higher for higher lens redshift: A6(/c,zl) cx 1/D(zl) 
(see Eq. ([2]))- Therefore, we find that even higher-z pop- 
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FIG. 15: Imprints of the local-type primordial non- 
Gaussianity in the galaxy-CMB lensing power spectrum, 
1(1 + 1)C; /(27r), for the same populations of lens galaxies as 
in Fig. [3] The solid, dashed, and dotted lines show /nl = 0, 
±50, and ±100, respectively. 



FIG. 16: Same as Fig. 1151 but for high-z lens galaxies with 
bi — 2 at zl = 2 (top-left), b\ — 2.5 at zl = 3 (top-right), 
b\ = 3 at zl = 4 (bottom-left), and i>i = 3.5 at zl — 5 
(bottom-right). 



ulations of galaxies give us a much better chance of de- 
tecting the effects of /nl- Figure [TBI shows C^ K (z L ) for 
/nl = ±50 and ±100 for populations of high-z galaxies: 
b\ = 2 at z L = 2 (top-left), b\ = 2.5 at zl = 3 (top- 
right), b\ = 3 at zl = 4 (bottom- left), and b\ = 3.5 at 
zl = 5 (bottom-right). The first one, a spectroscopic 
galaxy survey at Zl = 2 with 6^ = 2, is within reach by, 
e.g., the Hobby-Eberly Telescope Dark Energy Experi- 
ment (HETDEX) [11, |64|] . There we find, for / NL = ±50, 
~ 10% effect at I ~ 40, and a factor of two effect at I ~ 10 
(see Fig. ITT]) . The effects grow bigger at higher z: higher- 
z surveys at z > 3 can be done with, e.g., the concept 
of the Cosmic Inflation Probe (CIP) 5 . At z L = 4 and 5 
(with b\ = 3 and 4, respectively) we find ~ 10% effect at 
I ~ 100, a factor of two effect at I ~ 30, and even bigger 
effects at I < 30 (see Fig. Q7|). 



C. Covariance matrix of the galaxy-CMB lensing 

The covariance matrix of the galaxy-CMB lensing is 
given by .65j 

{C^ K C(t K ) - {C^){Cf, K ) 



(2/ + l)/ sky 



S W , (41) 



where JV" is the reconstruction noise from CMB given 
by [58| . The covariance matrix equation here is the same 
as Eq. ([M]) . except that now the shape noise of source 
galaxies is replaced by the reconstruction noise of CMB 
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FIG. 17: Same as Fig. 1111 but for high-z lens galaxies with 
bi — 2 at zl = 2 (top-left), b\ — 2.5 at zl = 3 (top-right), 
b\ = 3 at zl = 4 (bottom-left), and 6i = 3.5 at zl — 5 
(bottom-right). 



lensing. In what follows, we shall assume a "nearly per- 
fect" CMB experiment considered in Hu and Okamoto 
[58l|. w hose Gaussian random detector noise is modeled 
asm 



c 



_ ( T C MB \ 2 J(i+l) - 2 /81n2 



http:/ /www. cfa.harvard.edu/cip/ 



where the white noise level of detectors is At = 
Ap/V2 = 1 fiK arcmin, and the Full-Width-at-Half- 
Maximum (FWHM) of the beam is a = A'. With these 
detector parameters and the cosmological parameters of 
the "WMAP±BAO±SN ML" parameters in Table 1 of 
0, we find Ni ~ 6 x 10~ 8 sr" 1 on large scales, I < 100. 
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10 100 10 100 

Multipole, I Multipole, I 

FIG. 18: Angular power spectra of the galaxy-galaxy corre- 
lation, C; (thick dotted lines), the galaxy-convergence cross- 
correlation, C; (thick solid lines), and the convergence- 
convergence correlation, Cf (thick dashed lines) for the Gaus- 
sian initial condition (/nl = 0). The four panels show the 
same populations of galaxies and clusters of galaxies as in 
Fig. 1151 We also show the galaxy shot noise, 1/ul (thin dot- 
ted lines) as well as the lens reconstruction noise, Nj* (think 
dashed lines), for N L = 10 6 and ATf ~ 6 X 10" 8 sr" 1 (for 
multipoles much smaller than that corresponds to the beam 
size of 4'). We find l/n L < Cf and N? < C ; K for I < 100. 



Fig. [18] shows the galaxy-galaxy, galaxy-convergence, 
convergence-convergence angular power spectra for the 
Gaussian initial condition (/nl = 0). This figure is qual- 
itatively similar to Fig. [TSJ the galaxy-galaxy correlation 
is exactly the same, and the galaxy-convergence power 
spectrum is simply a scaled version of the corresponding 
curve in Fig. [12] The major difference comes from CP: 
as the CMB photons travel a longer path than photons 
from source galaxies, the convergence-convergence power 
spectrum is higher for the CMB lensing convergence. 

On large scales (I < 100), the covariance matrix is 
dominated by the cosmic variance terms: 1/ul -C C/ 1 
and Np -C Cp. The cross correlation coefficients are 
small, of order 10%: the maximum values are 0.12, 0.11, 
and 0.10 for Zl — 0.3, 0.5, and 0.8, respectively. There- 
fore, we can again use Eq. (|3"T)) for estimating the vari- 
ance, and find l m [ n (Eq. {38])) above which we can mea- 
sure the galaxy-convergence cross correlation with the 
signal-to-noise ratio greater than unity. For logarithmic 
bins of Al/l = 0.23, we find l min = 12.2, 13.5, and 15.8 
for zl = 0.3, 0.5, and 0.8, respectively. Comparing to the 
results in Sec IIII CI Z min is slightly bigger, as Cp (which 
contributes to the uncertainty) increases more rapidly 
than Cp K (the signal we are after) would as the source 
redshift increases from zs to z*. 

Fig. [TH] shows the expected 1-er uncertainties of the 
angular power spectrum of the galaxy-CMB lensing 
cross correlation, on top of the predicted Gaussian/non- 
Gaussian signals with five different values of non- 
Gaussianity parameters: /nl = 0, ±50, ±100. We also 
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FIG. 19: Same as Fig. 1151 but with 1-sigma uncertainty due 
to the shape noise of source galaxies (filled box, Eq. (129|) ) 
and full error budget (empty box, diagonal of Eq. (|27|ll in- 
cluding the cosmic variance. We use the multipole bins of 
size Al — 0.23Z. For uncertainty of CMB lensing reconstruc- 
tion, We assume the nearly-perfect reference experiment of 
Hu and Okamoto [k|: white detector noise At = Ap/\/2 = 
1 fiK arcmin, and FWHM of the beam a = 4'. 
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FIG. 20: Same as Fig. [18] but for the high redshift lens 
galaxies shown in Fig. 1161 For these populations (and with 
N L = 10 6 ), the shot noise is about the same as the galaxy 
power spectrum, i.e., C; ~ 1/ul. 



show the 1-a uncertainties without the cosmic variance 
due to the cosmic shear. Once again, it would be difficult 
to measure the effects of /nl from a single lens redshift, 
but combining many slices would help measure /nl from 
the galaxy-CMB lensing cross correlation. 

What about using even higher-z lens galaxies? As 
shown in Fig. 1201 for highcr-z populations (with z L = 
2 — 5) the galaxy-galaxy power spectra are about the 
same as the shot noise levels. This is true only for the 
assumed number of lenses, Nl = 10 6 (over the full sky), 
which is somewhat arbitrary. Increasing Nl will help 
reduce the noise, but only up to a factor of y/2. For 
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FIG. 21: Same as Fig. [19] but for the high redshift lens 
galaxies shown in Fig. 1161 



populations with Cf 



l/n-L, we can approximate the 



variance as 



Var (Cf" 



l/n L )C? 



2C/ l Cf 



(2Z + l)AZ/ sky (2; + l)A// sky " 



Thus, we find Cf K / Jv&y{C{ ik ) > 1 when 



l> 



r,V(Ai//)Aky 



(43) 



(44) 



The maximum cross-correlation coefficients are 0.091, 
0.084, 0.078, and 0.073 for z L = 2, 3, 4, and 5, re- 
spectively. The estimated / max is then 29 (zl = 2), 34 
(z L = 3), 38 {z L = 4) and 42 \z L = 5). 

In Fig. [21] we compare the expected 1-er uncertainties 
with the predicted signals from high-z lens galaxies with 
/nl = 0, ±50, and ±100. Comparing this result with 
that in Fig. [121 we conclude that higher- z lens popula- 
tions do provide a better chance of finding the effects 
of /nl than lower- z lenses, although we would still need 
to combine many lens redshift slices. In particular, us- 
ing highcr-z lenses, we can find non-Gaussian effects at 
higher and higher multipoles which are easier to measure; 
thus, high-z galaxies correlated with CMB lensing offers 
a yet another nice probe of the local-type primordial non- 
Gaussianity. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the galaxy-galaxy lens- 
ing and galaxy-CMB lensing cross-correlation functions. 
We have focused on large scales, typically larger than 
100 Mpc at the lens redshift. While current measure- 
ments have high signal-to-noise ratios on much smaller 
scales, we believe that future surveys will enable detec- 
tion of interesting physical effects in the large-scale, linear 
regime. 



We derive the full covariance matrix for galaxy-galaxy 
lensing, including the cosmic variance due to the clus- 
tering of lenses and to cosmic shear (Eq. [TH]). We use 
the linear bias model to provide the halo-mass and halo- 
halo correlations needed for this calculation. We present 
results for the covariance of the mean tangential shear 
measurement as a function of angular separations, as 
well as for the harmonic space halo-convergence cross- 
power spectrum. Our calculations show that the errors 
in AE(i?) are dominated by the cosmic variance term 
for R > 50 hr 1 Mpc (see Fig. [5]). Similarly, the errors 
in the halo-convergence cross power spectra, C^ K , are 
dominated by the cosmic variance term at I < 100 (see 

Fig. en 



For Gaussian initial conditions, we show that the bary- 
onic effects in the matter power spectrum (often called 
Baryon Acoustic Oscillations) produce a "shoulder" in 
the galaxy-galaxy lensing correlation (i.e., the mean tan- 
gential shears), AS(R), at R ~ 110 h' 1 Mpc (see Fig. [3]). 
This effect should be easy to measure from the next- 
generation lensing surveys by combining AS(i?) from 
multiple lens redshift slices. 

We consider the prospects of detecting primordial non- 
Gaussianity of the local-form, characterized by the /nl 
parameter. We have found that the scale-dependent 
bias from the local-form non-Gaussianity with /nl = 
±50 modifies AS(R) at the level of 10-20% at R ~ 
300 hr 1 Mpc (depending on b\ and z^; see Fig. [SJ (see 
Fig. [5]). The modification grows rapidly toward larger 
scales, in proportion to R 2 . High-z galaxies at, e.g., 
z > 2, cross-correlated with CMB can be used to find 
the effects of /nl in the galaxy-convergence power spec- 
trum, Cf K . While the effects are probably too small to 
see from a single lens redshift (see Fig. |2~T]) . many slices 
can be combined to beat down the cosmic variance er- 
rors. Exactly how many slices are necessary, or what is 
the optimal strategy to measure /nl from the galaxy- 
CMB lensing signal requires a more detailed study that 
incorporates the survey strategy for specific galaxy and 
lens surveys. 

We emphasize that, while the two-point statistics 
of shear fields are not sensitive to primordial non- 
Gaussianity, the two-point statistics correlating shear 
fields with density peaks (i.e., galaxies and clusters of 
galaxies) are sensitive due to the strong scale-dependence 
of halo bias on large scales. 

Finally, we note that one can also measure the effects 
of /nl on the halo power spectrum, Cj l . For example, C/ 1 
that would be measured from LSST can be used to probe 
/nl ~ 1 [13] ; thus, we would expect Cj 1 to be more pow- 
erful than the lens cross-correlation statistics we studied 
here. However, a combination of the two measurements 
would provide useful cross-checks, as galaxy clustering 
and galaxy-lensing correlations are affected by very dif- 
ferent systematics. 
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This completes the derivation of Eq. (flT|) . 
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APPENDIX A: DERIVATION OF THE MEAN 
TANGENTIAL SHEAR 



APPENDIX B: DERIVATION OF THE 
COVARIANCE MATRIX OF THE MEAN 
TANGENTIAL SHEAR 



One may write down the observed tangential shears at 
a given distance from a lens halo, 0, averaged over Nl 
lens halos as 



To compute the covariance matrix of the tangential 
shears (not yet azimuthally averaged), we first compute 



N L 

^2s D (n- hi) 



7t (n + 0), (Al) 



where Sd is the delta function, and i denotes the location 
of lens halos. Note that we have not azimuthally aver- 
aged the tangential shears yet. The ensemble average of 
7t l yields the number- weighted average of the tangential 
shear: 



(l' t l )(8) = ^-Jd 2 h(n L (h) lt (h + 0)), 



(A2) 



where nz,(n) is the surface number density of lens halos 
at a given location on the sky, n. Expanding it into the 
perturbation, n/,(n) = n/,[l + 5ft (n)], we obtain 



(lt){0) 



1 



/sky 



d 2 h 

47T 



(S h (h) lt (h + 6)), 



(A3) 



where / s k y = N^/ (Ann^) is a fraction of sky covered by 
the observation. From statistical isotropy of the universe, 
(5ft(n)7 t (n + 0)) does not depend on n, and thus the 
integral over n simply gives 47r/ s k y . Expanding 5h and 
7t in Fourier space, we obtain 



(lt){0) 



d 2 \ d 2 Y 

(2^(2^ 
d 2 \ 



e il.4 e fl -(n+fl) cos [ 2 (<f,-p)}(6 h (l) K (l')) 
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C^cos[2(^-^)] ( 
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where we have used (S h (Y)K(Y)) = (2tt) 2 C/ ik ( 5 d (1 + !')■ 
Finally, we take the azimuthal average of (^) (6) to find 



x (S D (h - hi)S D (ti - n 3 ) 7t (n + 0) lt (ti + 0')) 
' ' d 2 h I d 2 h' 



N 2 

x [6 D (h - h')(n L (h)j t (h + 6)j t (h' + 0')) 

+(n L (h)n L (h') Jt (h + 6) lt {ti + 6'))] . (Bl) 

Here, the first term in the square bracket correlates 
two 7t's measured relative to the same lens halo (1-halo 
term), and the second correlates two 74 's relative to two 
lens halos (2-halo term). Again expanding til into the 
perturbation, riL(n) = Ul[1 + 5ft (n)], we obtain 



( 7t (0) 7t (0')) 

T~i!r I ^(it(h + e) lt (h + 9>)) 

/sky N L J 4tt 
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d 2 h f d 2 h' 
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47T 



[< 7t (n + e) 7t (n' + e')> 



+{5 h (h)5 h (h') 7t (h + e) lt {ti + e'))\ .(B2) 

Here, we assume that Sh and 7 t obey Gaussian statistics, 
i.e., ((5ft 7 t7t) = 0. This approximation is justified even 
in the presence of primordial non-Gaussianity, as non- 
Gaussianity is weak, and this approximation only affects 
the size of errorbars. Let us evaluate each term. With 7t 
expanded in Fourier space, the first term (1-halo term) 
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becomes 

1 1 



mean tangential shear: 
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(cf K ) 2 + cf + 




where <7 7 is the r.m.s. shape noise (dimensionless) , and 
ng is the surface density of source (background) galaxies 
that are available for the shear measurement at a given 
location. By azimuthally averaging 7t , we find 



This completes the derivation of Eq. ([T 
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APPENDIX C: ON THE ACCURACY OF 
LIMBER'S APPROXIMATION 
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Here, C ; K is the angular power spectrum of 

As for the second term (2-halo term), the first of the 
second term vanishes, as J ePn 7 t(n + 9) = 0. The re- 
maining non- vanishing term gives 
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Throughout this paper we have repeatedly used Lim- 
ber's approximation in order to relate the angular cor- 
relation function to the corresponding three dimensional 
power spectrum. In general, Limber's approximation is 
(B4) known to be accurate only for small angular scales, and 
only for the quantities which are integrated over a broad 
range of redshift. 

However, the situations we have considered in this pa- 
per sometimes violate both of the conditions above: 1) 
We correlate the convergence field with galaxies within a 
very thin redshift slice, and 2) the non-Gaussianity signal 
we study in this paper appears only on very large scales. 

Then, how accurate is Limber's approximation in this 
case? In this Appendix, we shall study in detail the va- 
lidity and limitation of Limber's approximation, by com- 
paring the main results of the paper to the result of exact 
calculations. 

Consider a quantity Xj(n), which is projected on the 
sky. Here, n is the unit vector pointing toward a given 
n (e-e') direction on the sky. This quantity is related to the three 
dimensional quantity Sj(r; z) by a projection kernel Wi(z) 
as 



(B5) 



Here, C/ 1 is the angular power spectrum of 5^(1). By 
azimuthally averaging 7t in the above equation, we find 



Xi(n) = dzW l (z)s i [d A {z)n; z] 
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where we have used the identity 
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Throughout this Appendix, we use d A (z) to denote 

d A (0;z). 

Fourier transforming Sj(r), one obtains 

Si[d A (z)h; z)} 
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Collecting both the 1-halo and 2-halo terms, we finally 
obtain the covariance matrix of the azimuthally-averaged 



In the third line, we have used Rayleigh's formula: 
e *W = 4^^ i J J; (fcr)y ; ; i (k)T im (n). 



l,m 
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By using Eq. (|C2[) , we rewrite Eq. (|C1[) as 
Xi(h) = 4?r^V / dzW l {z) 



l,m 

d 3 k 
J2nf 



Si (k, 2 )i i [fcd A (z)]Y^(k)F im (n](.C3) 



Therefore, the coefficients of the spherical harmonics de- 
composition of Xi(xi), a**, becomes 

d 3 k 



a Z = 4 ™ / dzWi(z) 



(27T) 



3 s 4 (k,z)j z [fcd A (z)]^(k). 

(C4) 

We calculate the angular power spectrum, Cf ,x ' , by tak- 
ing an ensemble average of (»fm a Zm*) as 
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P s ^(k;z,z r ) 
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(2tt) 3 

xi,[fcdA(^)]ii[fcdA(«')]^(k)^(k). 
where we have used the definition of the power spectrum: 

( Si (k, z)s*(k', z)) = (2nf6(k - k')P s > s > (k; z, z'). 

Now, by assuming statistical isotropy of the universe, 
we write P SiS ' (k; z, z') = P SiS i (fc; z, z'), and do the angu- 
lar integration of k by using the orthonormality condition 
of spherical harmonics: 



dkY* m (k)Y4(k) = 1. 
We then obtain the angular power spectrum given by 
Cf iXi = J dzWi{z) j dz'Wj{z') 

k 2 dkP s ^ (k; z, z% [kd A {z)] 3l [kd A (z')} 



is a sharply peaked function, we cannot use Limber's ap- 
proximation given below. This is the reason why we have 
used the exact result for C/\ 

In order to get the expression for Limber's approxima- 
tion, we assume that P SiS i (k) is a slowly-varying function 
of k. Then, by using the identity 
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we approximate the k integral of Eq. (|C6p as 
- / k 2 dkP s ^ (k)ji(kr)ji(kr') 
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By using this approximation, we finally get 
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which is the result known as Limber's approximation. 

One important application of Limber's approximation 
is the statistics involving weak gravitational lensing. The 
lensing kernel for the convergence field, W K (z), can be 
calculated by integrating the lens equation: 



W K (z) 



Po 



Y, c (z;z s )H(z) ' 
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where T, c (z; zg) is the critical surface density defined in 
Eq. (fT6|) . The exact result for the galaxy-convergence 
angular cross power spectrum is 



C^(z L ) = - 



dz- 



Po 



7r J T, c (z;z s )H(z) 
x J dkk 2 P hm (k,z L , z)ji[kd A {z L )]ji[kd A {z)] 1 

(C13) 



This is the exact relation. 

What determines the form of Wi(z)l For a projected 
galaxy distribution projected on the sky, this kernel is 
simply a normalized galaxy distribution function in red- 
shift space. In this paper, we consider the delta function- 
like distribution, i.e., 



(C6) and the exact result for the convergence-convergence an- 
gular power spectrum is 



dz 
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x dktff > P m {k,z,z%[kd A {zM[kdA{z')]. (C14) 



W g (z) = S D (z-z L ). 



First, we compare the exact convergence-convergence 

(C7) 

angular power spectrum to Limber's approximation. 

TT • t — i r^rrry .,, , , ,, „ ,. , , . Fig. [22l shows that Limber's approximation works very 

Using Eq. (IC6I) with the delta function kernel above „ r „ r , , .,, , , 

. , , ° well tor all tour source redsnilts we study m tne paper: 

yields Eq. K3BJ): ' ' 



C[ l = - / dkk 2 P g (k, z L )jf [kd A {z L )\ . (C8) 



z s = 0.65, 1.19, 2.25, and 1089.0. For I > 10, the er- 
ror caused by Limber's approximation is always much 
smaller than 1%. 

Then, we compare the galaxy-convergence cross angu- 
Again, this is still the exact result. As the form of lar power spectra. Fig. |2"31 and Fig. [2^1 show the compar- 
W g (z) we have considered here (i.e., a delta function) ison between the exact galaxy-convergence cross power 
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FIG. 22: Top: Convergence-convergence angular power 
spectrum from two different methods: the exact calculation 
(Eq. lC14l symbols) and Limber's approximation (Eg. 1351 solid 
lines). Bottom: Fractional differences between Limber's ap- 
proximation and the exact integration. Symbols are the same 
as the top panel. Grey symbols show the absolute values of 
negative values. 




: z =0.3, z =0.65 

L S 

: z L =0.5, z =1.19 
: z =0.8, z =2.25 



FIG. 24: Same as Fig. 1221 but for the non-Gaussian cor- 
rection (i.e., the term proportional to Ab(k)) to the galaxy- 
convergence cross angular power spectrum. We show the cor- 
rections with /nl = 1 and fei = 2. 



10~ J 
lO" 4 

-? 10" 5 

CM 

> 10- 6 

- 10" 7 

+ 

^ 10" 8 

10- 9 

7 io° 
1 io-' 

S"_ 10" 2 
o 

> 10" 3 
1 10" 4 



CMB lensing 


Non-biosed (b ] 


= 1), Goussion (f =0) ^^a^^m 






■^^^£5^ ■ z =0-3 






■^^^^^^ "■ \=0.5 






„ : z =0.8 _ 






»: z|=2.0 






□ : z =3.0 - 






x : z =4.0 






. : z =5.0 


H 1 


" l 8 ■ a . 
« . * a 

* f a 


°* • * * * " ■ * : 



10 

multipole, I 



100 



10 

multipole, I 



100 



FIG. 25: Same as Fig. [23] but for the galaxy-CMB lensing. 



FIG. 23: Same as Fig. I22I but for the galaxy-convergence 
cross angular power spectrum with /nl = and bi = 1. 



spectrum (Eq. IC13| symbols) and their Limber approxi- 
mation (Eq. [321 solid lines) for three galaxy-galaxy lens- 
ing cases we study in Sec. [TTJ (zl,Zs) = (0.3,0.65), 
(0.5,1.19), and (0.8,2.25). 

For the Gaussian term (Fig. I23p . Limber's approxima- 
tion is accurate at I > 10 with the errors less than 1%. 
On the other hand, Limber's approximation to the non- 
Gaussian correction term (Fig. [24]) has a sizable error, at 
the level of 10%, at I ~ 10. The error goes down to the 
1% level only at I ~ 100. One needs to keep this in mind 
when comparing Limber's approximation with observa- 



tions. We find that Limber's approximation underpre- 
dicts the Gaussian term at I < 20, while it overpredicts 
the non-Gaussian corrections at all multipoles. 

The story is basically the same for the galaxy-CMB 
lensing cross power spectrum. Fig. [25] (Gaussian term) 
and Fig. [21] (non-Gaussian correction) show the compar- 
ison between the exact galaxy-convergence cross power 
spectrum (Eq. IC13[ solid lines) and their Limber approx- 
imation (Eq. [32] dashed lines) for seven lens redshifts we 
study in SecHV] Zl — 0.3, 0.5, 0.8, 2, 3, 4, and 5. Again, 
for small scales, I > 10, Limber's approximation works 
better than 1% for the Gaussian term, while it overpre- 
dicts the non-Gaussian correction at the level of 10% at 
I ~ 10 and 1% at I ~ 100. 

What about the effect on the mean tangential shear, 
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FIG. 26: Same as Fig. [24] but for the galaxy-CMB lensing. 
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FIG. 27: Top: Same as Fig. [3] but also showing the exact 
result (Eq. IC13I thick lines) on top of the result from Lim- 
ber's approximation (Eq. I32I thin lines). Bottom: Fractional 
difference of Limber's approximation relative to the exact re- 
sult. 

AS (J?)? Fig. [27] compares the Gaussian term of AS(i?) 
from the exact integration and that from Limber's ap- 
proximation. On the top panel of Fig. 12 71 we show 
the baryonic feature computed with Limber's approxi- 
mation (thin lines, the same as those in Fig. [3]) as well as 
that computed with the exact integration (thick lines). 
They are indistinguishable by eyes. The bottom panel 
shows the fractional differences between the two. We 
find that Limber's approximation is better than 0.5% for 
R < 180 hT x Mpc; thus, the baryonic feature in AS is 
not an artifact caused by Limber's approximation. 

However, Limber's approximation becomes worse and 
worse as we go to larger R. Fig. [21] shows AS on large 
scales. For the lens redshifts that we have studied here, 



FIG. 28: Same as Fig. [27] but for larger R. Thick lines are 
the results of the exact integration, while the thin lines are 
Limber's approximation. The Limber approximation over- 
predicts AE(_R) for large R, but the error is at most 5% for 
R < 500 h~ Mpc. The error is the largest for the lowest zl, 
as a physical separation R at a lower redshift corresponds to 
a larger angular separation on the sky. 
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FIG. 29: Fractional differences in the non- Gaussian cor- 
rection terms, AE n G, from Limber's approximation and the 
exact integration. Using Limber's approximation, we over- 
predict the non-Gaussian correction by ~ 20% at R = 
300 h' 1 Mpc for z L = 0.3. 



the error is at most 5% for R < 500 h^ 1 Mpc, and the 
error is the largest for the lowest zl, as a given R at a 
lower redshift corresponds to a larger angular separation 
on the sky. 

While Limber's approximation undcrpredicts the 
Gaussian term on large scales, it overpredicts the non- 
Gaussian correction terms. Fig. fJH shows the fractional 
differences of the non-Gaussian correction terms, AS n G, 
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FIG. 30: Same as Fig. [4] but with the exact integration 
instead of Limber's approximation. 



between Limber's approximation and the exact calcula- 
tion as a function of separation R for three lens redshifts: 
zl = 0.3, 0.5, and 0.8. This figure shows that the error 
caused by Limber's approximation can be substantial on 
AE nG . 

As Limber's approximation to AE(i?) can be quite in- 
accurate on very large scales, we show the exact calcu- 
lations of AE(i?) in Fig. [30] (Limber's approximation is 
given in Fig. 01) 

Finally, we note that the definition of the tangential 
shear we have used (Eq. |4|) is valid only on the flat sky 
(as noted in the footnote there), and thus the prediction 
for AS on very large scales probably needs to be revisited 
with the exact definition of the tangential shears on the 
full sky using the spin-2 harmonics. This is beyond of 
the scope of our paper. 
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